The man, the plant, and the insect: shooting host specificity determinants in Serratia marcescens pangenome.

Introduction: Serratia marcescens is most commonly known as an opportunistic pathogen causing nosocomial infections. It, however, was shown to infect a wide range of hosts apart from vertebrates such as insects or plants as well, being either pathogenic or growth-promoting for the latter. Despite being extensively studied in terms of virulence mechanisms during human infections, there has been little evidence of which factors determine S. marcescens host specificity. On that account, we analyzed S. marcescens pangenome to reveal possible specificity factors. Methods: We selected 73 high-quality genome assemblies of complete level and reconstructed the respective pangenome and reference phylogeny based on core genes alignment. To find an optimal pipeline, we tested current pangenomic tools and obtained several phylogenetic inferences. The pangenome was rich in its accessory component and was considered open according to the Heaps' law. We then applied the pangenome-wide associating method (pan-GWAS) and predicted positively associated gene clusters attributed to three host groups, namely, humans, insects, and plants. Results: According to the results, significant factors relating to human infections included transcriptional regulators, lipoproteins, ABC transporters, and membrane proteins. Host preference toward insects, in its turn, was associated with diverse enzymes, such as hydrolases, isochorismatase, and N-acetyltransferase with the latter possibly exerting a neurotoxic effect. Finally, plant infection may be conducted through type VI secretion systems and modulation of plant cell wall synthesis. Interestingly, factors associated with plants also included putative growth-promoting proteins like enzymes performing xenobiotic degradation and releasing ammonium irons. We also identified overrepresented functional annotations within the sets of specificity factors and found that their functional characteristics fell into separate clusters, thus, implying that host adaptation is represented by diverse functional pathways. Finally, we found that mobile genetic elements bore specificity determinants. In particular, prophages were mainly associated with factors related to humans, while genetic islands-with insects and plants, respectively. Discussion: In summary, functional enrichments coupled with pangenomic inferences allowed us to hypothesize that the respective host preference is carried out through distinct molecular mechanisms of virulence. To the best of our knowledge, the presented research is the first to identify specific genomic features of S. marcescens assemblies isolated from different hosts at the pangenomic level.

For Someone, You Are the Whole World: Host-Specificity of Salmonella enterica.

Salmonella enterica is a bacterial pathogen known to cause gastrointestinal infections in diverse hosts, including humans and animals. Despite extensive knowledge of virulence mechanisms, understanding the factors driving host specificity remains limited. In this study, we performed a comprehensive pangenome-wide analysis of S. enterica to identify potential loci determining preference towards certain hosts. We used a dataset of high-quality genome assemblies grouped into 300 reference clusters with a special focus on four host groups: humans, pigs, cattle, and birds. The reconstructed pangenome was shown to be open and enriched with the accessory component implying high genetic diversity. Notably, phylogenetic inferences did not correspond to the distribution of affected hosts, as large compact phylogenetic groups were absent. By performing a pangenome-wide association study, we identified potential host specificity determinants. These included multiple genes encoding proteins involved in distinct infection stages, e.g., secretion systems, surface structures, transporters, transcription regulators, etc. We also identified antibiotic resistance loci in host-adapted strains. Functional annotation corroborated the results obtained with significant enrichments related to stress response, antibiotic resistance, ion transport, and surface or extracellular localization. We suggested categorizing the revealed specificity factors into three main groups: pathogenesis, resistance to antibiotics, and propagation of mobile genetic elements (MGEs).

Recombination in Bacterial Genomes: Evolutionary Trends.

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Salmonella-Based Biorodenticides: Past Applications and Current Contradictions.

The idea of using pathogens to control pests has existed since the end of the 19th century. Enterobacteria from the genus Salmonella, discovered at that time, are the causative agents of many serious diseases in mammals often leading to death. Mostly, the strains of Salmonella are able to infect a wide spectrum of hosts belonging to vertebrates, but some of them show host restriction. Several strains of these bacteria have been used as biorodenticides due to the host restriction until they were banned in many countries in the second part of the 20th century. The main reason for the ban was their potential pathogenicity for some domestic animals and poultry and the outbreaks of gastroenteritis in humans. Since that time, a lot of data regarding the host specificity and host restriction of different strains of Salmonella have been accumulated, and the complexity of the molecular mechanisms affecting it has been uncovered. In this review, we summarize the data regarding the history of studying and application of Salmonella-based rodenticides, discuss molecular systems controlling the specificity of Salmonella interactions within its multicellular hosts at different stages of infection, and attempt to reconstruct the network of genes and their allelic variants which might affect the host-restriction mechanisms.

Current Methods for Recombination Detection in Bacteria.

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

Identification of Genetic Risk Factors of Severe COVID-19 Using Extensive Phenotypic Data: A Proof-of-Concept Study in a Cohort of Russian Patients.

The COVID-19 pandemic has drawn the attention of many researchers to the interaction between pathogen and host genomes. Over the last two years, numerous studies have been conducted to identify the genetic risk factors that predict COVID-19 severity and outcome. However, such an analysis might be complicated in cohorts of limited size and/or in case of limited breadth of genome coverage. In this work, we tried to circumvent these challenges by searching for candidate genes and genetic variants associated with a variety of quantitative and binary traits in a cohort of 840 COVID-19 patients from Russia. While we found no gene- or pathway-level associations with the disease severity and outcome, we discovered eleven independent candidate loci associated with quantitative traits in COVID-19 patients. Out of these, the most significant associations correspond to rs1651553 in MYH14p = 1.4 × 10-7), rs11243705 in SETX (p = 8.2 × 10-6), and rs16885 in ATXN1 (p = 1.3 × 10-5). One of the identified variants, rs33985936 in SCN11A, was successfully replicated in an independent study, and three of the variants were found to be associated with blood-related quantitative traits according to the UK Biobank data (rs33985936 in SCN11A, rs16885 in ATXN1, and rs4747194 in CDH23). Moreover, we show that a risk score based on these variants can predict the severity and outcome of hospitalization in our cohort of patients. Given these findings, we believe that our work may serve as proof-of-concept study demonstrating the utility of quantitative traits and extensive phenotyping for identification of genetic risk factors of severe COVID-19.

Temporal Control of Seed Development in Dicots: Molecular Bases, Ecological Impact and Possible Evolutionary Ramifications.

In flowering plants, seeds serve as organs of both propagation and dispersal. The developing seed passes through several consecutive stages, following a conserved general outline. The overall time needed for a seed to develop, however, may vary both within and between plant species, and these temporal developmental properties remain poorly understood. In the present paper, we summarize the existing data for seed development alterations in dicot plants. For genetic mutations, the reported cases were grouped in respect of the key processes distorted in the mutant specimens. Similar phenotypes arising from the environmental influence, either biotic or abiotic, were also considered. Based on these data, we suggest several general trends of timing alterations and how respective mechanisms might add to the ecological plasticity of the families considered. We also propose that the developmental timing alterations may be perceived as an evolutionary substrate for heterochronic events. Given the current lack of plausible models describing timing control in plant seeds, the presented suggestions might provide certain insights for future studies in this field.

Dissecting the Environmental Consequences of Bacillus thuringiensis Application for Natural Ecosystems.

Bacillus thuringiensis (Bt), a natural pathogen of different invertebrates, primarily insects, is widely used as a biological control agent. While Bt-based preparations are claimed to be safe for non-target organisms due to the immense host specificity of the bacterium, the growing evidence witnesses the distant consequences of their application for natural communities. For instance, upon introduction to soil habitats, Bt strains can affect indigenous microorganisms, such as bacteria and fungi, and further establish complex relationships with local plants, ranging from a mostly beneficial demeanor, to pathogenesis-like plant colonization. By exerting a direct effect on target insects, Bt can indirectly affect other organisms in the food chain. Furthermore, they can also exert an off-target activity on various soil and terrestrial invertebrates, and the frequent acquisition of virulence factors unrelated to major insecticidal toxins can extend the Bt host range to vertebrates, including humans. Even in the absence of direct detrimental effects, the exposure to Bt treatment may affect non-target organisms by reducing prey base and its nutritional value, resulting in delayed alleviation of their viability. The immense phenotypic plasticity of Bt strains, coupled with the complexity of ecological relationships they can engage in, indicates that further assessment of future Bt-based pesticides' safety should consider multiple levels of ecosystem organization and extend to a wide variety of their inhabitants.

The Distribution of Several Genomic Virulence Determinants Does Not Corroborate the Established Serotyping Classification of Bacillus thuringiensis.

Bacillus thuringiensis, commonly referred to as Bt, is an object of the lasting interest of microbiologists due to its highly effective insecticidal properties, which make Bt a prominent source of biologicals. To categorize the exuberance of Bt strains discovered, serotyping assays are utilized in which flagellin serves as a primary seroreactive molecule. Despite its convenience, this approach is not indicative of Bt strains' phenotypes, neither it reflects actual phylogenetic relationships within the species. In this respect, comparative genomic and proteomic techniques appear more informative, but their use in Bt strain classification remains limited. In the present work, we used a bottom-up proteomic approach based on fluorescent two-dimensional difference gel electrophoresis (2D-DIGE) coupled with liquid chromatography/tandem mass spectrometry(LC-MS/MS) protein identification to assess which stage of Bt culture, vegetative or spore, would be more informative for strain characterization. To this end, the proteomic differences for the israelensis-attributed strains were assessed to compare sporulating cultures of the virulent derivative to the avirulent one as well as to the vegetative stage virulent bacteria. Using the same approach, virulent spores of the israelensis strain were also compared to the spores of strains belonging to two other major Bt serovars, namely darmstadiensis and thuringiensis. The identified proteins were analyzed regarding the presence of the respective genes in the 104 Bt genome assemblies available at open access with serovar attributions specified. Of 21 proteins identified, 15 were found to be encoded in all the present assemblies at 67% identity threshold, including several virulence factors. Notable, individual phylogenies of these core genes conferred neither the serotyping nor the flagellin-based phylogeny but corroborated the reconstruction based on phylogenomics approaches in terms of tree topology similarity. In its turn, the distribution of accessory protein genes was not confined to the existing serovars. The obtained results indicate that neither gene presence nor the core gene sequence may serve as distinctive bases for the serovar attribution, undermining the notion that the serotyping system reflects strains' phenotypic or genetic similarity. We also provide a set of loci, which fit in with the phylogenomics data plausibly and thus may serve for draft phylogeny estimation of the novel strains.

Analysis of the Spectrum of ACE2 Variation Suggests a Possible Influence of Rare and Common Variants on Susceptibility to COVID-19 and Severity of Outcome.

OBJECTIVES: In March 2020, the World Health Organization declared that an infectious respiratory disease caused by a new severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2, causing coronavirus disease 2019 (COVID-19)] became a pandemic. In our study, we have analyzed a large publicly available dataset, the Genome Aggregation Database (gnomAD), as well as a cohort of 37 Russian patients with COVID-19 to assess the influence of different classes of genetic variants in the angiotensin-converting enzyme-2 (ACE2) gene on the susceptibility to COVID-19 and the severity of disease outcome. RESULTS: We demonstrate that the European populations slightly differ in alternative allele frequencies at the 2,754 variant sites in ACE2 identified in the gnomAD database. We find that the Southern European population has a lower frequency of missense variants and slightly higher frequency of regulatory variants. However, we found no statistical support for the significance of these differences. We also show that the Russian population is similar to other European populations when comparing the frequencies of the ACE2 variants. Evaluation of the effect of various classes of ACE2 variants on COVID-19 outcome in a cohort of Russian patients showed that common missense and regulatory variants do not explain the differences in disease severity. At the same time, we find several rare ACE2 variants (including rs146598386, rs73195521, rs755766792, and others) that are likely to affect the outcome of COVID-19. Our results demonstrate that the spectrum of genetic variants in ACE2 may partially explain the differences in severity of the COVID-19 outcome.

No More Tears: Mining Sequencing Data for Novel Bt Cry Toxins with CryProcessor.

Bacillus thuringiensis (Bt) is a natural pathogen of insects and some other groups of invertebrates that produces three-domain Cry (3d-Cry) toxins, which are highly host-specific pesticidal proteins. These proteins represent the most commonly used bioinsecticides in the world and are used for commercial purposes on the market of insecticides, being convergent with the paradigm of sustainable growth and ecological development. Emerging resistance to known toxins in pests stresses the need to expand the list of known toxins to broaden the horizons of insecticidal approaches. For this purpose, we have elaborated a fast and user-friendly tool called CryProcessor, which allows productive and precise mining of 3d-Cry toxins. The only existing tool for mining Cry toxins, called a BtToxin_scanner, has significant limitations such as limited query size, lack of accuracy and an outdated database. In order to find a proper solution to these problems, we have developed a robust pipeline, capable of precise 3d-Cry toxin mining. The unique feature of the pipeline is the ability to search for Cry toxins sequences directly on assembly graphs, providing an opportunity to analyze raw sequencing data and overcoming the problem of fragmented assemblies. Moreover, CryProcessor is able to predict precisely the domain layout in arbitrary sequences, allowing the retrieval of sequences of definite domains beyond the bounds of a limited number of toxins presented in CryGetter. Our algorithm has shown efficiency in all its work modes and outperformed its analogues on large amounts of data. Here, we describe its main features and provide information on its benchmarking against existing analogues. CryProcessor is a novel, fast, convenient, open source (https://github.com/lab7arriam/cry_processor), platform-independent, and precise instrument with a console version and elaborated web interface (https://lab7.arriam.ru/tools/cry_processor). Its major merits could make it possible to carry out massive screening for novel 3d-Cry toxins and obtain sequences of specific domains for further comprehensive in silico experiments in constructing artificial toxins.

Phenome-wide functional dissection of pleiotropic effects highlights key molecular pathways for human complex traits.

Over the recent decades, genome-wide association studies (GWAS) have dramatically changed the understanding of human genetics. A recent genetic data release by UK Biobank (UKB) has allowed many researchers worldwide to have comprehensive look into the genetic architecture of thousands of human phenotypes. In this study, we used GWAS summary statistics derived from the UKB cohort to investigate functional mechanisms of pleiotropic effects across the human phenome. We find that highly pleiotropic variants often correspond to broadly expressed genes with ubiquitous functions, such as matrisome components and cell growth regulators; and tend to colocalize with tissue-shared eQTLs. At the same time, signaling pathway components are more prevalent among highly pleiotropic genes compared to regulatory proteins such as transcription factors. Our results suggest that protein-level pleiotropy mediated by ubiquitously expressed genes is the most prevalent mechanism of pleiotropic genetic effects across the human phenome.

Reproductive History of a Woman With 8p and 18p Genetic Imbalance and Minor Phenotypic Abnormalities.

We report on the phenotype and the reproductive history of an adult female patient with an unbalanced karyotype: 8p23 and 18p11.3 terminal deletions and 8p22 duplication. The indication for karyotyping of the 28-year-old patient was a structural rearrangement in her miscarriage specimen: 45,ХХ,der(8;18)t(8;18)(p23;p11.3). Unexpectedly, the patient had the same karyotype with only one normal chromosome 8, one normal chromosome 18, and a derivative chromosome, which was a product of chromosomes 8 and 18 fusion with loss of their short arm terminal regions. Fluorescence in situ hybridization revealed that derivative chromosome was a pseudodicentric with an active centromere of chromosome 8. Array comparative genomic hybridization confirmed 8p and 18p terminal deletions and additionally revealed 8p22 duplication with a total of 43 OMIM annotated genes being affected by the rearrangement. The patient had minor facial and cranial dysmorphia and no pronounced physical or mental abnormalities. She was socially normal, had higher education and had been married since the age of 26 years. Considering genetic counseling, the patient had decided to conceive the next pregnancy through in vitro fertilization (IVF) with preimplantation genetic testing for structural chromosomal aberrations (PGT-SR). She underwent four IVF/PGT-SR cycles with a total of 25 oocytes obtained and a total of 10 embryos analyzed. Only one embryo was balanced regarding chromosomes 8 and 18, while the others were unbalanced and demonstrated different combinations of the normal chromosomes 8 and 18 and the derivative chromosome. The balanced embryo was transferred, but the pregnancy was not registered. After four unsuccessful IVF/PGT-SR cycles, the patient conceived naturally. Non-invasive prenatal testing showed additional chromosome 18. The prenatal cytogenetic analysis of chorionic villi revealed an abnormal karyotype: 46,ХХ,der(8;18)t(8;18)(p23;p11.3)mat,+18. The pregnancy was terminated for medical reasons. The patient has a strong intention to conceive a karyotypically normal fetus. However, genetic counseling regarding this issue is highly challenging. Taking into account a very low chance of balanced gametes, emotional stress caused by numerous unsuccessful attempts to conceive a balanced embryo and increasing age of the patient, an IVF cycle with a donor oocyte should probably be considered.

Whole-exome sequencing provides insights into monogenic disease prevalence in Northwest Russia.

BACKGROUND: Allele frequency data from large exome and genome aggregation projects such as the Genome Aggregation Database (gnomAD) are of ultimate importance to the interpretation of medical resequencing data. However, allele frequencies might significantly differ in poorly studied populations that are underrepresented in large-scale projects, such as the Russian population. METHODS: In this work, we leveraged our access to a large dataset of 694 exome samples to analyze genetic variation in the Northwest Russia. We compared the spectrum of genetic variants to the dbSNP build 151, and made estimates of ClinVar-based autosomal recessive (AR) disease allele prevalence as compared to gnomAD r. 2.1. RESULTS: An estimated 9.3% of discovered variants were not present in dbSNP. We report statistically significant overrepresentation of pathogenic variants for several Mendelian disorders, including phenylketonuria (PAH, rs5030858), Wilson's disease (ATP7B, rs76151636), factor VII deficiency (F7, rs36209567), kyphoscoliosis type of Ehlers-Danlos syndrome (FKBP14, rs542489955), and several other recessive pathologies. We also make primary estimates of monogenic disease incidence in the population, with retinal dystrophy, cystic fibrosis, and phenylketonuria being the most frequent AR pathologies. CONCLUSION: Our observations demonstrate the utility of population-specific allele frequency data to the diagnosis of monogenic disorders using high-throughput technologies.